version 12.0
clear
capture log close
set more off

*     ****************************************************************  *
*     ****************************************************************  *
*       Filename:       PSED_data_feature_analysis.do                   *
*       Date:           06/01/2014                                      *
*       Author:         Martin Ottmann                                  *
*       Input File(s):  PSED_data_feature_analysis.dta                  *
*     ****************************************************************  *
*     ****************************************************************  *


// only PROMISES
logit pcmEND pcm1 pcm1_squared pcm1_cubed ///
      incomp durationDY intDY unpko multireb ///
	  ppsPROM mpsPROM epsPROM tpsPROM, ///
	  robust cluster(pseddyID)
estimates store feature1

// only PROMISED IMPLEMENTATION
logit pcmEND pcm1 pcm1_squared pcm1_cubed ///
      incomp durationDY intDY unpko multireb ///
	  ppsPROM_IMP mpsPROM_IMP epsPROM_IMP tpsPROM_IMP, ///
	  robust cluster(pseddyID)
estimates store feature2

// PROMISES and PROMISED IMPLEMENTATION
logit pcmEND pcm1 pcm1_squared pcm1_cubed ///
      incomp durationDY intDY unpko multireb ///
	  ppsPROM mpsPROM epsPROM tpsPROM ///
	  ppsPROM_IMP mpsPROM_IMP epsPROM_IMP tpsPROM_IMP, ///
	  robust cluster(pseddyID)	  
estimates store feature3

estout       ///
   feature1  ///
   feature2  ///
   feature3, ///
   cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) starlevels(* 0.10 ** 0.05 *** 0.01) stats(aic bic ll N)
